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Abstract 

We study the thermal phase transition in colour SU(N C = 3) Quantum Chromodynamics (QCD) with a variable number 
of fermions in the fundamental representation by using lattice Monte-Carlo simulations. We collect the (pseudo) 
critical couplings /3 L C for Nf = (0, 4, 6, and 8), and we investigate the pre-conformal dynamics associated with the 
infra-red fixed point in terms of the Nf dependence of the transition temperature. We propose three independent 
estimates of the number of flavour Nf where the conformal phase would emerge, which give consistent results within 
the largish errors. We consider lines of fixed N, in the space of (Nf, bare lattice coupling), and locate the vanishing 
of the step scaling function for Nf ~ 11.1 ± 1.6. We define a typical interaction strength gj at the scale of critical 
temperature T c and we find that g T c meets the zero temperature critical couplings estimated by the two-loop Schwinger 
Dyson equation or the IRFP coupling in the four-loop beta-function at Nf ~ 12.5 + 0.7. Further, we study the Nf 
dependences of T c /M where M is a UV Nf independent reference scale determined by utilising the coupling at the 
scale of the lattice spacing. Then, T c /M turns out to be a decreasing function of Nf and the vanishing T c /M indicates 
O ' the emergence of the conformal window at N* f ~ 10.4 + 1 .2. 
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1. Introduction 

The analysis of the phases of strong interactions present many fascinating challenges - mechanisms of the confine- 
ment, a different realisation of the chiral symmetry, nature of the symmetric phase, the emergence of the conformality, 
and so on. All these topics are under active scrutiny both theoretically and experimentally QJJ] . While strong interac- 
tions spontaneously break a chiral symmetry in the ordinary QCD at zero temperature, the chiral symmetry is realised 
either at high temperatures - in the so called quark-gluon plasma (QGP) phase - and at a large number of flavours 
Nf > N^ (even at zero temperature) J2, S HI . In the latter case, the theory is expected to become not only chirally 
but also conformally invariant. This is due to the emergence of an infra-red fixed point (IRFP) for Nf > Nf at a cou- 
pling which is not strong enough to break the chiral symmetry. Both physics intuition and phenomenological analysis 
based on functional renormalisation group ffl and finite temperature holographic QCD J3l indicate that the conformal 
phase of cold, many flavour QCD and the high temperature chirally symmetric phase are continuously connected. 
In particular, the onset of the conformal window coincides with the vanishing of the transition temperature, and the 
conformal window appears as a zero temperature limit of a possibly strongly interacting QGP. 

The analysis of the finite temperature phase transition is a well established line of research within the lattice 
community, and our approach will be completely conventional here. According to the Pisarski-Wilczek scenario (H], 
the most likely possibility for Nf > 3 is a first order chiral transition in the chiral limit, turning into a crossover 
above a critical mass endpoint, and/or on lattices which are not large enough. We will identify such crossover with 
confidence for a number of flavours ranging from four to eight, and we will complement these results with those of 
the deconfinement transition in the quenched model. Then, we study the approach to the conformal phase in the light 
of the chiral phase transition at finite temperature with variable number of flavours. 

One problem of this approach is the setting of a common scale among theories which are essentially different. 
We will propose two alternative possibilities to handle this problem, one evolving from our previous work Jgt], and 
the other from a recent analysis 111 Oil . Interestingly, this latter approach analyses the dependence of the confinement 
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parameters on the matter content, and proposes microscopic mechanisms for confinement motivated by such ty 
dependence. Further, we will argue that even results in the bare lattice parameters can be used directly to locate the 
critical number of flavours, thus generalising to finite temperature the Miransky-Yanawaki phase diagram, Ref. yfl. 
A second zero of the two-loop beta-function of a non-Abelian gauge theory implies, at least perturbatively, the 



appearance of IRFP conformal symmetry 111 ll 11211 . In colour SU(3) gauge theory with ty massless fundamental 



fermions, the second zero appears at ty- > 8.05, before the loss of asymptotic freedom (LAF) at ty/ AF = 16.5. Ana- 
lytic studies of the conformal transition of strong interactions have produced a variety of predictions for the conformal 
threshold: the Schwinger-Dyson approach with rainbow resummations J2, S 01 or the functional renormalisation 
grou p method J6J suggest the onset of conformal window around N*. ~ 12. An all-order perturbative beta-function 



1 13| inspired by the NSVZ beta-function of SQCD yj] leads to a bound N* f > 8.25. Instanton studies at large ty 11511 
claimed a qualitative change of behavior at ty = 6. ty has also been estimated for different fermion representations 



116[. 

The sub-critical region, when Nf gets closer and closer to Np is interesting per-se: the question is whether the 
chiral dynamics there shows any difference with the standard QCD dynamics. Significant differences with respect 
to the QCD dynamics might offer a basis to model builders interested in beyond-the-standard-model theories. The 
recent discovery of a 125 GeV boson at the LHC poses the question as to whether there are light composite scalars 
which might be identified with such boson, in alternative to a standard model Higgs boson. Pre-conformal dynamics 
might well help these studies J5[[l3l. In our study, such pre-conformal dynamics could manifest itself either with a 
clear observation of a separation of scales, or with a manifestation of a critical behaviour when approaching ty. One 
possibility is to observe the Miransky-Yamawaki essential singularity [3]. Alternatively, in an FRG approach Is], the 
pseudocritical line is almost linear with Nf for small Nf, and displays a singular behaviour when approaching ty, 
which could be the only observable effects, beyond Miransky scaling. A jumpin g sc enario in which the change from 
a QCD dynamics to the conformal window is abrupt is also a distinct possibility ll 811 . 

Clearly, as in any system undergoing a phase transition, the nature and extent and of critical window are a purely 
dynamical question whose answer cannot be guessed a priori. Since the underlying dynamics is completely non- 
perturbative, lattice calculations are the onl y to ol to perform an ab-initio, rigorous study of these phenomena, and 



many lattice studies have recently appeared IU9I1 



This paper is one step of our ongoing program 12011 - 112811 which aims at elucidating the phase diagram of QCD 



with fundamental fermions on the lattice, and in the continuum. Further studies either with fundamental fermions \2S 
- lf39ll or other representations J40I l4ll I42I I43I Hi have contributed to our current understanding of this challenging 



field. However, only a subset of these studies has addressed issues related with pre-conformal dynamics 14511 - 114811 . 
i2, 36, 44 1 which are the main theme of this paper. The direct inspection of theories at fixed Nf is often inconclusive, 
especially close to the expected threshold ty*.. Also because of this, we feel it is an useful approach to try to observe 
directly the approach to conformality by monitoring the evolution of the pre-conformal results as a function of Nf. 

In this paper, we investigate the thermal chiral phase transition for Nf = 0,4,6,8 colour SU(N C = 3) QCD by 
using lattice QCD Monte Carlo simulations with sta gger ed fermions. Here, Nf — 6 and 8 is expected to be in the 
important regime as sug gested by the results in Refs. 1 15, 29ll . We combine our findings with those of our early work 



forty = 6 and 8 P, |28ll . 

This paper grows out of our early study {si and extends it in several ways: We have accumulated more statistics 
and added more parameters, and we present here an extended set of simulations and details. We develop a new scale 
setting procedure, so that we can more confidently measure the critical temperature on a common reference scale 
among different number of flavour. Furthermore, we present new estimates of the critical number of flavours ty: 
Partly motivated by the recent work [[loll , we introduce a typical interaction strength g£ at the critical temperature 
based on our lattice results, and compare it with a four-loop IRFP coupling (#irfp,4i) 14911 and a critical coupling 
C?sd) estimated by using a two-loop Schwinger-Dyson equation 115011 . Further, we introduce and discuss the finite 
temperature version of the Miransky-Yamawaki phase diagram, and propose a strategy to locate the critical number 
of flavour motivated by the properties of the lattice step scaling function in the vicinit y of the IRFP i34ll . Some of the 
new results presented here have been anticipated in a recent proceeding, and talks Il20ll2 111 . 

This paper is organized as follows: In the next section, we will explain the simulation setups. In Section [3] we 
show our results for the chiral crossover at finite T, for each ty, and then, we collect the critical lattice couplings 
associated with the chiral crossovers at ty = 0, 4, 6, 8. In Section |4] we investigate the asymptotic scaling of our 
critical couplings at each ty . In Section [5] we investigate the ty dependences of the chiral crossovers, and estimate 
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Figure 1: Left: The Monte-Carlo trajectories of the chiral condensate (PBP) obtained by using the lattice volume 24 3 X 12 just before the chiral 
crossover /Sl = 5.45 - 5.50 at Nf = 6. Right: The jackknife errors as a function of a bin-size for the trajectories shown in the left panel. 



the lower edge of the conformal window N^. Finally in Section [6] we provide concluding remarks. The Appendix 
is devoted to the summary tables of the simulation parameters and the numerical results obtained by analysing the 
simulation outputs. 



2. Simulations' setup 

We investigate finite temperature QCD with different number of flavour Nf = (0, 4, 6, 8) by utilising the publicly 
available MILC code [51]. The temperature T is related to the inverse of the lattice temporal extension. 



T = . (1) 

a(fiL)-N, 

and control it by varying /?l at fixed N t . The number of lattice points in the spacial directions N s is chosen such that 
the aspect ratio N s /N, > 2 in all our runs. For each Nf, we use a single bare fermion mass ma = 0.02. The simulation 
parameters used in this study are summarised in | Appendix A 



2.1. Action and algorithm 

The setup for the action explained below is the same as the one used for Nf — 8 in Ref. lEsll up to the number 
of flavour. We use an improved version of the staggered action, the Asqtad action ll52Tl . with a one-loop Symanzik 



1531 154|] and tadpole [55!!] improved gauge action, 



-^-Tr log M[am, [/,«<,] + £ MglW - U Q ] , (2) 



where g^ is the lattice bare coupling, and /J, are defined as 

x /10 jS„(l -0.4805or s ) B„ \ 
Q3 P ,f3 r ,/3 pg ) = - - PpK 2 ^ -^f0.03325a. s (3) 

a s = -4 log — ^— , ito = (U c > 1/4 ■ (4) 
s 3.0684 x " 

The plaquette coupling fi p = lO/g^ = Pl is a simulation input. The M[am, U, uq] in Eq. (f2]i denotes the matrix for a 
single flavour Asqtad fermion with bare lattice mass am, and Uq represents the trace of the ordered product of link 
variables along C,, for the lxl plaquettes (i = p), the 1x2 and 2x 1 rectangles (i = r), and the lxlxl parallelograms 
(i = pg), respectively - all divided by the number of colours. The tadpole factor uq is determined by performing zero 
temperature simulations on the 12 4 lattice (the second column of Table I A. 101 - IA. 14l i, and used as an input for finite 
temperature simulations. 



3 



To generate configurations with mass degenerate dynamical flavours, we have used the rational hybrid Monte 
Carlo algorithm (RHMC) J5^], which allows to simulate an arbitrary number of flavours Nf through varying the 
number of pseudo-fermions. The quenched (Nf = 0) system has been realised by using massive bare fermion mass 
ma = 1 .0 in the four flavour system. The six flavour system has been realised by using two pseudo-fermions in the 
Rational approximation with a quarter root technique, Nf =4-2-3/4 = 6. Then, we have assumed the rooting does 
not affect the results within the accuracy of our simulation. For the other number of flavour (Nf = 0, 4, 8), we do not 
use the rooting. 

We have adjusted the microcanonical step length St and the step length of a single trajectory At = 20 x St to 
realise 75 - 80 percent Metropolis acceptances. Details are reported in the fourth column of Table IA. 10l - IA. 1 4l For 
each parameter set, we have collected a number of trajectories ranging from a one thousand to ten thousands - the 
latter closer to the chiral crossover regime. 

2.2. Observables 

The focus of this paper is the analysis of the chiral phase transition The fundamental observable is then the order 
parameter for chiral symmetry, the chiral condensate: 

fl3<#> = 4^ Tr[M " ] ) ' (5) 

where N s (N t ) represents the number of lattice sites in the spatial (temporal) direction. We have measured a 3 (if/if/) by 
using a stochastic estimator with 20 repetitions. We have also measured connected and disconnected chiral suscepti- 
bilities, 

^Xccmn = --^-(Tt\(MM)-'\) , 

4NsN t L J 



m 

2 / 

a *di s " 



\6N]N, 



(Tr [m -1 ] 2 ) - (Tr [m -1 ]) 2 



(6) 



Here we have conveniently written the chiral condensate and its susceptibilities in terms of traces of (products of) the 
staggered fermion matrix M. We note that the MILC convention for the chiral condensate gives the twice of Eq. 10, 
as will be indicated several times in the following sections for results. We have measured the susceptibilities a 2 ^ C onn 
and fl 2 ^disc separately. 

The disconnected chiral susceptibility is a non-local quantity which can be estimated from the variance of the 
bulk behaviour of the chiral condensate. Since we have used the stochastic estimator for the chiral condensate mea- 
surements, the variance would automatically include part of the connected contributions through random sources 



multiplying themselves. Following Bernard et al. 15711 . we take account this effect in our estimate for the disconnected 
part a 2 x&sc by considering the only off-diagonal elements of the covariance matrix for the random sources. 

The measurements of a 3 (i(/4r) and a 2 ^f C onn,diss allow us to construct two physically relevant quantities: the scalar 
and pseudo-scalar susceptibilities, 

= _ ^ ™ 

Xo~ — ^ — Xconn ' ^disc ? \ ' ) 

om 

Xn = - (8) 

m 

Their associated cumulant 

Rn=~, (9) 

Xn 



is a probe of the chiral symmetry [28, 58]- This is owing to the fact that^o- and;^ are related through Ward identities to 
the spacetime volume integral of the scalar (cr) and pseudoscalar (n) propagators. In the chiral limit, the susceptibility 
ratio R„ should be one in chirally symmetric regime due to the degeneracy of the chiral partners, while it should be 
zero in the spontaneously broken phase. Even including a finite bare fermion mass, R„ still has a strong signal for 
the chiral transition or crossover. In particular, R„ ~ 1.0 in the chirally symmetric regime holds true till the chiral 
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Table 1: Summary of the (pseudo) critical lattice couplings y3 L c for the theories with Nf = 0, 4, 6, 8, am = 0.02 and varying N t = 4, 6, 8, 12. 
The entries with * are the update for our previous results (jj. The entries with f have been quoted from our previous studies on Nf = 8 (2^1 . 



N f \N t 


4 


6 


8 


12 





7.35 + 0.05 


7.97* ± 0.07 


8.26 + 0.06 




4 


5.65 + 0.05 


6.00* ± 0.05 


6.15 + 0.15 




6 


4.675* + 0.05 


5.025* +0.05 


5.20* ± 0.05 


5.55* + 0.1 


8 




4.1 125"'" + 0.0125 


4.275 ± 0.05 


4.34 1 " + 0.04 



condensate is dominated by the linear mass term contribution. It turns out that R K allows the identification of a pseudo 
critical coupling /3 L C associated with the chiral crossover, which , in the cases we have studied, coincides in the error 
with the pseudocritical coupling determined from the maximum of the chiral susceptibility. 
In the gauge sector, we measure the Polyakov loop, 



^^z^nM' (io) 



where tr e denotes the trace in colour space, and U$ <t *. is the temporal link variable. From the variance of L, we also 
evaluate the susceptibility for the Polyakov loops. 
2.3. Statistics, and error analysis 

In the vicinity of the chiral crossover, we have a long auto-correlation time, and thermalization checks require 
extra care. Here we explain our analyses by using a typical example: The left panel of Fig. [T]displays the evolution of 
the chiral condensate on a the lattice volume 24 3 x 12 just before the chiral crossover p^ = 5.45 and 5.50 at Nf = 6, 
one of the most time-consuming examples in our simulations. (In order to shorten the simulation time, we started 
the evolution from thermalised configurations obtained at /3l < 5.45.) We have computed the ensemble averages by 
using the last 2500 (2000) trajectories at /?l = 5.45 (5.50) and we have confirmed that they are consistent with those 
obtained by using last 2000 (1500) trajectories. We have then used the latter trajectories to evaluate the average. In 
the cases we are considering, the correspond to the data found in the right hand side of the vertical green (dashed) 
lines in the left panel of Fig. Q] 

We divide the obtained data set into several bins and utilise the jackknife method in order to take account of the 
auto-correlation effect into the error estimate. As a bin-size Sbin becomes larger, the jackknife error increases (the 
right panel of Fig. [O, which is due to the decrease of the effective number of (uncorrelated) data (« ave / Sbin, "ave = the 
number of trajectories to calculate the average). For a sufficiently large Sbin, the jackknife errors at/?L = 5.45 and 5.50 
level off, giving a reliable error estimate. 

Here is the result obtained from the above procedures: 

Ph «traj "ave «bin 2a 3 (ifr^) 

5.45 3100 2000 400 0.0622(2), 

5.50 4100 1500 375 0.0570(5), 

We have performed the analyses explained here for all the various p^, Nf, and the lattice volumes. The results are 



summarised in Appendix A 



3. Results on the lattice thermal transition 

In this section, we show our simulation results on the chiral and deconfinement crossover for the different number 
of flavours Nf. 

We have used a common bare fermion mass ma = 0.02 for all simulations at finite Nf. According to the Pisarski- 
Wilczek scenario 18[], the most likely possibility for Nf > 3 is a first order chiral transition in the chiral limit. In- 
troducing a bare fermion mass, the first order phase transition will eventually turn into a crossover for masses larger 

5 



than some critical mass. Since the chiral condensate looks smooth in our results, we are most likely above the critical 
endpoint in all the cases we have studied, and we use the terminology of "chiral crossover" in the following. 

The finite bare mass ma = 0.02 might have a different physical relevance at each Nf, as well as for different 
bare coupling for a fixed Nf. It remains then to be see how our results would change in the chiral limit, and we 
hope to come back to this point in a future study. Since we have noted that at strong coupling and small masses the 
improvement term in the Action might be resposible for the spurious phases 123l 12211 observed also in Ref. 132l 12511 
we might favor an unimproved action for this study. 

Before entering into details, let us summarise our main results, i.e. the critical lattice couplings /? L C associated with 
the chiral crossover in Table Q] For Nf = (4, 6, 8) we have observed that the peak position of the chiral susceptibility 
a 2 x<T, whenever clearly defined, coincides within the errors with the inflection point of R K defined in Eq. (|9), as well as 
with the inflection point of the chiral condensate and that of the Polyakov loop. This indicates that the crossover region 
is rather narrow, as different indicators give consistent pseudocritical points. We then quote the common pseudocrical 
coupling, with a conservative error estimate. For the quenched (Nf = 0) case, we have extracted the pseudocritical 
coupling from the deconfinement crossover by extracting the peak position of the Polyakov loop susceptibility. 

In the following subsections, we present these results in detail, starting from Nj — 6 and 8 in the first two 
subsections, and continuing with the Nf = 4 and Nf = 0. The reader who is not interested in these technical details is 
advised to skip the rest of this Section and proceed directly to the next one. 



3.1. Chiral crossover at Nf — 6 

We show the Nf = 6 results for a fixed bare fermion mass ma = 0.02. In Figs. [2] the chiral condensate a 3 {i[/i[/) 
(PBP, red O) the real-part of Polyakov loop L (Re[PLOOP], blue □), the chiral susceptibility (Xa-, red +), and the 
chiral susceptibility ratio (R„, blue x) are displayed as a function of a lattice coupling /?l = 10/g£. The first, second, 
third, and fourth lines in the figure show the results obtained by using temporal extensions N, = 4, 6, 8, and 12, 
respectively. We shall now extract the critical lattice couplings /? L C associated with the thermal chiral crossover from 
these results. 

As shown in the left panel of the first line in Fig. the largest decrease of chiral condensates (PBP, red O) (as 
well as a drastic increase of the real part of Polyakov loops (RefPLOOP], blue □)) is found in between /?l = 4.65 and 
4.70. Thus, we expect the chiral crossover in this region. As shown in the right panel of the first line in Fig. 01, the 
chiral susceptibility a 2 Xcr (red +) has a clear peak at /3l = 4.65. The Gaussian fit for the susceptibilities in the range 
[4.4, 4.9] leads to a maximum at a slightly larger /?l (red dashed line). Further, the susceptibility ratio R„ (blue x) has 
an inflection point around /?l = 4.65 - 4.70. For larger the increasing rate of R„ significantly reduces, and gets 
to almost unity. Thus, all observables consistently indicate the pseudo-critical coupling to be y6 L c = 4.675 + 0.05 for 
(Nf, N t ) = (6, 4). The error is determined to include the next to neighbour data and the maximum of the Gaussian fit. 

The second line in Fig. [2]i displays the results for N, = 6. As shown in the left panel, the largest decrease of chiral 
condensates (PBP, red O) (as well as a drastic increase of the real part of Polyakov loops (RefPLOOP], blue □)) is 
found in between /?l = 5.00 and 5.05, and we expect the chiral crossover in this region. As shown in the right panel, 
the chiral susceptibility a 2 Xcr (red +) has a peak at /?l = 5.05, and the Gaussian fit for the susceptibilities in whole 
range of /?l has a maximum at a slightly smaller /3l = 5.0 (red dashed line). The susceptibility ratio R„ (blue x) has 
an inflection point around /3l = 5.00 - 5.05, and then, it goes into the plateau domain. All observables consistently 
indicate the pseudo-critical coupling to be /? L C = 5.025 + 0.05 for (Nf, N t ) = (6, 6). The error is determined to include 
both jSl = 5.0 and 5.05 enough. 

The third line in Fig. |2]i shows the results for N, = 8. As indicated by the left panel, the chiral condensates as 
well as the Polyakov loops look smooth at almost everywhere, and it is difficult to locate the crossover point from 
them. As shown in the right panel, the chiral susceptibility a 2 x<r (red +) has a peak at /?l = 5.2, and the Gaussian fit 
for the susceptibilities in whole range of /?l has a maximum at a slightly smaller /?l = 5.17 (red dashed line). The 
susceptibility ratio R n (blue x) exhibits the largest variation between /3l = 5.15 and 5.2, after which the increasing 
rate of R„ reduced and eventually evolves into almost unity. From the peak position of a 2 Xa-, we estimate the critical 
coupling to be around = 5.20 + 0.05 for (Nf,N t ) = (6, 8). The error is determined to include the next neighbour 
data, the maximum of the Gaussian fit for^ tr , and the R„ inflection point. 

Finally, we analyse the results for N, = 12, the largest temporal extension in our Nf — 6 simulations. The /?l 
dependence of chiral condensates are found to be particularly smooth for whole range of /?l = 4.7 - 5.7. Note that 



6 



0.7 
0.6 hi 
0.5 
0.4 
0.3 
0.2 
0.1 



@ © 


®© 5 








2xPBP h©h 
RerPI.OOPl hFH 










© 



C 
















) 


















© 




Izl 














© 


Id 

3 












1- 
















f 


© 


s 


I l-LI 11=] 1*1 1XJ 










* 



3.5 3.75 4 4.25 4.5 4.75 5 5.25 5.5 

Pl 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 












2xPBP h©h 
4xRerPLOOPl >» 



























© 




a 










c 






a 














3 
















V 














1 3 E 






Itl "J 


KB 









4.25 4.5 4.75 

Pl 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



© 












2xPBP 


l-0H 












x KelPLUUFI 






































6 






















© 










UJ 








c- 
























V 


















B ©0 










: R : 


3 L 










5 , <=> 















4.6 4.7 4.8 4.9 





8e-02 r 




7e-02 - 




7e-02 - 


- 


6e-02 - 


a 




E 


6e-02 - 








5e-02 - 




5e-02 - 




4e-02 - 




4e-02 



©© 



5 5.1 5.2 5.3 5.4 5.5 

Pl 

4e-03 
3e-03 
3e-03 
2e-03 
2e-03 
le-03 



O 
O 
J 

3S 



5e-04 
0e+00 
-5e-04 



5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 



3 
2.5 

2 
1.5 

1 

0.5 


















WlorXo -- 


H 






























( 4 














/ 




\ - 


, *L 


x+ ! 






* ■ 


yf 


' > 

X 


: \ y 








< 


v--<f X > 


1 
< 











4 4.2 4.4 4.6 4.8 5 5.2 5.4 

Pl 



2 - 



tit lor x„ 

Xa H-h 



/f 



^ 



xx V 



xxx 



.XX 



3.6 3.8 4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 

Pl 



1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 


1.2 
1.1 
1 

0.9 
0.8 











u 
















;. 








V, 












X 
















> 


x 

X 

< 




lit for x„ — - 


) 


< * X 






X H-h " 



4.6 4.8 



5 5.2 5.4 

Pl 



5.6 



i i 






































.T. X 




< 










> 


c 














X 5 












~i 


I 

E 


C > 


i 
















5 


la H-h 
















1 



5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 
Pl 



Figure 2: The Nt = 6 results for a fixed bare fermion mass ma = 0.02. The first, second, third, and fourth lines show the results obtained by using 
temporal extensions N, = 4, 6, 8, and 12, respectively. In each line, the left panel shows the chiral condensate in lattice unit (PBP, red O) and 
the real-part of Polyakov loops (RefPLOOP], blue □), and the right panel displays the chiral susceptibility (Xa-, red +) and the chiral susceptibility 
ratio (R„, blue x), as a function of jSl- For N t = 4, the Gaussian fit for the chiral susceptibility has been performed in the range [4.4, 4.9] to capture 
the peak structure. 
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Figure 3: The iV/ = 8 results obtained by using 24 3 X 8 lattice volume with ma = 0.02: The chiral condensate in lattice unit (PBP, red Q in left 
panel), the real-part of Polyakov loop (Re[PLOOP], blue □ in left panel), and the chiral susceptibility (Xa-, red + in right panel), and the chiral 
susceptibility ratio (R K , blue X in right panel) are shown as a function of 



in this case, the aspect ratio (N s /N t ) is only two, and larger volumes would be required to reach a comparable clarity 
in the signal. As shown in the left panel of final line in Fig. the onset for the Polyakov loop at /?l = 5.525 is still 
appreciable (blue □). We here notice that the increase of Polyakov loops so far has been found just before the chiral 
crossover in the case of smaller temporal extensions N, =4-8, though the Polyakov loop itself is not associated 
with the chiral dynamics. Based on such an experience, we assume that the chiral crossover at N, = 12 is in the 
vicinity of the onset of the Polyakov loop, and carefully investigate the corresponding region /3l = 5.35 - 5.60. The 
chiral condensates do not have any clear signal (red O in the left panel). As shown in the right panel, the chiral 
susceptibility a 1 Xa (red +) has a small peak-like structure at (3^ = 5.575, and a bump-like structure at /?l = 5.50. 
The chiral susceptibility ratio R„ (blue x) has the largest increase between /?l = 5.55 and 5.575, and tends to be fiat 
in /?l > 5.575. Thus, the critical lattice coupling would be in the range 5.50 < jS L c < 5.575. Here, we employ a 
conservative estimate /? L C =5.55 + 0.1, which sufficiently covers the whole candidate range. 

Our y8 L c collection at Nf = 6 is found in the third line of Table [TJ As will be shown in the next subsection, the 
N, dependent nature of /? L C at Nf — 6 (a thermal scaling) is associated with the uniqueness of the physical critical 
temperature, indicating the chiral (non-conformal) dynamics at Nf = 6. 



3.2. Chiral crossover at Nf — 8 

We shall now investigate Nf = 8 system. In our previous paper II 2 811 . we have studied the chiral phase transition 
at Nf — 8 by using two lattice temporal extensions: N, — 6 and 12. One of the main results was that the chiral phase 
transition at Nf = 8 still showed a thermal scaling property, which indicated the existence of a typical scale associated 
with the chiral dynamics rather than the conformality. We here add additional data computed at N, = 8, and confirm 
the thermal scaling at Nf = 8, for this largish mass. 

The left panel of Fig. [3] shows ensemble averages of chiral condensates a 3 (ifnfr) (PBP, red OX the real-part of 
Polyakov loop L (Re[PLOOP], blue □) as a function of /3l- We observe the largest decrease of the chiral condensate 
between /?l = 4.25 and 4.30, while the real part of the Polyakov loop stars growing around /3l = 4.25. Although the 
error is huge, the chiral susceptibility ratio R n seems to have a larger increase between /?l = 4.25 and 4.30. The large 
error of R„ at /3l = 4.25 comes from the very long correlation in the Monte-Carlo trajectories, and hints at criticality. 
All observations consistently indicates the critical coupling to be (3£ = 4.275 + 0.05. Combining with N, — 6 and 12 
data 12811 . we summarise the critical coupling y6 L c at Nf — 8 in the final line of Table [TJ 



Here we should put some caveats on the Nf — 8 results: First, we have not observed a clear peak-like structure 
in the chiral susceptibility Xv- We should probably study larger volumes, with similar aspect ratio as those studied 
in other cases. The location of the pseudocritical point might change. Second, even in the strong coupling region 
P < j8 L c , R„ shows relatively large value ~ 0.7 - 0.8 We should go to even stronger coupling and smaller masses before 
observing a clear mass gap. And third, the N, dependence of /? L C shows a large deviation from the two-loop asymptotic 
scaling law as will be shown in the next subsection. This is not surprising, given that the couplings explored for 
Nf = 8 are larger than in other cases. This could imply something which cannot be captured at two-loop, for example, 
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the pre-conformal dynamics. Apparently, these caveats call for more detailed and quantitative lattice studies with a 
larger lattice size and a smaller bare fermion mass before drawing definite conclusions on Nf = 8. We note a recent 
study claiming the conformality emerges for Nf — 8 for small enough quark masses IjDl . 

3.3. Chiral crossover at Nf — 4 

Now we extend our analyses into the Nf = 4 system. In Figs.|H we show the chiral condensate a 3 (i/fi/f) (PBP, red 
O) the real-part of Polyakov loop L (PLOOP, blue □), the chiral susceptibility (x, T , red +), and the chiral susceptibility 
ratio (R n , blue x) are displayed as a function of a lattice coupling /?l = 10/g£. The first, second, and third lines in the 
figure show the results obtained by using temporal extensions N, = 4, 6, and 8, respectively. 

As shown in the left panel of the first line in Fig. [4] the largest decrease of chiral condensates (PBP, red O) (as 
well as a drastic increase of the real part of Polyakov loops (Re f PLOOP], blue □)) is found in between /3l - 5.60 
and 5.65, and, we expect the chiral crossover in this region. As shown in the right panel of the first line in Fig. |4j 
the chiral susceptibility a 2 Xcr (red +) gets to a maximum at /?l = 5.65. The Gaussian fit for the susceptibilities in the 
range [5.4, 5.8] leads to a maximum at a slightly larger /?l (red dashed line). Further, the susceptibility ratio R„ (blue 
x) has an inflection point around /?l = 5.60 - 5.65. For larger the increasing rate of R n significantly reduces, and 
eventually evolves into almost unity. Thus, all observables consistently indicates the pseudo-critical coupling to be 
/? L C = 5.65 + 0.05 for (Nf, N t ) = (4, 4). The error is determined to include the next to neighbour data and the maximum 
of the Gaussian fit. 

The second line in Fig. |4]i displays the results for N, = 6. As shown in the left panel, the chiral condensates 
(PBP, red O) are found to be smooth, and it is difficult to locate the chiral crossover. The real part of Polyakov loops 
(RefPLOOP], blue □) starts increasing around /?l = 5.95. Based on our previous experiences, the chiral crossover 
could be around this region. As shown in the right panel, the chiral susceptibility a 2 Xa (red +) has a maximum at 
/?l = 6.00, and the Gaussian fit for the susceptibilities in the range of [5.8,6.2] has a maximum at /?l = 6.00 (red 
dashed line). From the maximum position of the chiral susceptibilities, we estimate the pseudo-critical coupling to be 
j8 L ° = 6.00+0.05 for (Nf,N) = (4,6). The error is determined to include the next-to-neighbour data. The susceptibility 
ratio R n (blue x) has a significant increase in fi^ > 5.9, and then, it goes into the plateau domain at /?l = 6.1. This 
behaviour would be consistent to the above estimate /? L C = 6.00. 

The third line in Fig.|4]shows the results for N, = 8. As indicated by the left panel, the chiral condensates as well 
as the Polyakov loops look smooth at almost everywhere, and it is difficult to locate the crossover point from them. 
As shown in the right panel, the chiral susceptibility a 2 Xcr (red +) has a peak at /?l = 5.95. Indeed the susceptibility 
ratio R n (blue x) also shows a bump structure around /?l = 5.95, and implies some kinds of an instability of the 
system. However, the value of R n at /3l = 5.95 turns out to be at most 0.4, which indicates a large remaining of 
the chiral symmetry breaking. In turn, R„ keeps increasing till it approximately reaches to unity at fi^ = 6.30, and 
thus, the first peak at fi^ = 5.95 would not well capture the position of the chiral crossover. We here postpone the 
precise determination of the chiral crossover, and just provide a rough estimate of the pseudo critical coupling: The 
susceptibility ratio R„ has a large increasing rate in the range [6.0, 6.3]. We adopt the intermediate value as the pseudo 
critical coupling with the error covering whole range of [6.0, 6.3],y6 L c = 6.15 + 0.15. This also includes the maximum 
of the Gaussian fit for the chiral susceptibility /?l = 6.04. 

The second line of Table [TJprovides a summary of y6 L c for Nf = 4. 

3.4. Deconfinetnent at Nf — 

In this subsection, we estimate the critical lattice coupling j3£ for deconfinement in the quenched (Nf = 0) system. 
We note that both deconfinement and chiral transitions are associated with the thermal phase transition from the 
hadronic phase to the non-Abelian plasma phase with a drastic increase of the pressure (degrees of freedom). Then, 
our interest is the probe of the system with various Nf in light of such a thermal phase transition. In this sense, we 
regard the deconfinement crossover at Nf — as a continuation to the chiral crossover at finite Nf. In our setup, 
this connection is made explicit by the fact that we are realising a quenched system by use of a large mass in the 
four flavour system. It should be noted, however, that our result for the (pre-)conformal dynamics does not crucially 
depend on the quenched data. 

In Figs. [5] the thermalised ensemble averages of the absolute of Polyakov loop (\L\, blue □), its susceptibility (x\l\, 
red +), are displayed as a function of a lattice coupling j8l = 10/g?. The first, second, and third lines in the figure 
show the results obtained by using temporal extensions N, = 4, 6, and 8, respectively. 
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Figure 4: The Nf = 4 results for a fixed bare fermion mass ma = 0.02. The first, second, and third lines show the results obtained by using temporal 
extensions N t = 4, 6, and 8, respectively. In each line, the left panel shows the chiral condensate in lattice unit (PBP, red O) and the real-part of 
Polyakov loops (RefPLOOP], blue □), and the right panel displays the chiral susceptibility (Xa-, red +) and the chiral susceptibility ratio (R^, blue 
x), as a function of /3 L - For N, = 4 and 6, the Gaussian fits for the chiral susceptibilities have been performed in the range [5.4, 5.8] and [5.8, 6.2], 
respectively. 
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Figure 5: The quenched results. The first, second, and third lines show the results obtained by using temporal extensions N t = 4, 6, and 8, 
respectively. In each line, the left panel shows the absolute values of Polyakov loops (|L|, □), and the right panel displays the susceptibility 
calculated from the variance of \L\ (Xm> symbol +). The Gaussian fits are particularly bad here, however the identification of the maximum is clear. 
For Nt = 8, the fit has been performed excluding the data in /3 L > 8.8, because as shown in the left panel, the drastic increase of \L\ is found for a 
much smaller /? L - 
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As shown in the left panel of the first line in Fig. [5] the largest increase of the absolute value of Polyakov loops 
(\L\, blue □)) is found in between /?l = 7.30 and 7.35, and we expect the deconfinement crossover in this region. 
As shown in the right panel, the susceptibility for \L\ (red +) has a clear peak at /3l = 7.35, hence we estimate the 
pseudo-critical coupling to be = 7.35 + 0.05 for (N/,N t ) = (0,4). The error is determined to include the next to 
neighbour data. 

The second line in Fig. [5] displays the results for N, = 6. As shown in the left panel, the largest increase of 
the absolute value of Polyakov loops (|L|, blue □)) is found in between /?l = 7.80 and 7.90, and we expect the 
deconfinement crossover in this region. The maximum of the susceptibility evaluated from \L\ (red +) is observed at 
Nf = 7.9. The large error indicates the long correlation time of the Monte-Carlo trajectories. The Gaussian fit for the 
susceptibility has a maximum at /?l = 7.97. From this, we estimate the pseudo-critical coupling to be /? L C = 7.97 + 0.07 
for (Nf, Nt) = (0, 6). The error is determined to include the next-to-next-to-neighbour data from the maximum point. 

The third line in Fig.[5]represents the results for N, = 8. As shown in the left panel, the absolute value of Polyakov 
loops (\L\, blue □)) start increasing around /?l = 8.15, and we expect the deconfinement crossover in this region. The 
maximum of the susceptibility evaluated from \L\ (red +) is observed at Nf = 8.25. Again, we find the large error in 
the vicinity of the maximum, indicating the long correlation time of the Monte-Carlo trajectories. The Gaussian fit for 
the susceptibility in the range [7.8, 8.6] has a maximum at /?l = 8.26. Note that the fit range sufficiently covors whole 
region of the drastic increases of \L\ shown in the left panel. We adopt the maximum of the Gaussian fit as a critical 
coupling: = 8.26 + 0.06 for (Nf,N,) = (0,6). The error is determined to include the next-to-next-to-neighbour 
data from the maximum point. 

The first line of Table [Uprovides a summary of /? L C for Nf = 0. The N, dependent nature of y6 L c reflects the thermal 
nature of the crossover. 



4. Asymptotic scaling analyses for chiral phase transition 

In this section, we investigate the asymptotic scaling of the pseudocritical temperatures T c 

where /? L C have been computed in the previous section, and discuss the connection to the continuum physics. In 
the first subsection 14. II we introduce the normalised critical temperature 7\./Al/e JHoll from our y6 L c , where Al (Ae) 
represents the lattice (E-scheme) Lambda-parameter defined in the two-loop perturbation theory with or without a 
renormalisation group inspired improvement l6lll . Then in the subsections l4.2l the asymptotic scaling will be assessed 
by studying the (in)dependence of Tc/Al/e at each Nf. 

4.1. Normalised critical temperature 7" c /Al/e 
We here consider the two-loop beta function 

p(g) = -(bog 3 +b lg 5 ), (12) 

1 ( HC 2 [G] 4T[F]N f \ 
b °=(4n?[— 3— J' (13) 

* = ^pf^ - + , d4) 

with(C2[G], C2[F], T[F]) = (N c , (N^ — l)/(2N c ), 1/2). The coupling g is now assumed to be the lattice bare coupling 
gh = or the E-scheme renormalised coupling gE = -\/3(l - (P)(gi)), where {P)(g]_) is a plaquette measured 

in the zero temperature simulation with the input /Jl = 10/g L - In particular, the E-scheme critical coupling g E c is 
obtained by using our critical coupling g L c = -sJlO/fi^: g E c = -J3(l - (P)(g^)). If the one-loop perturbation theory 
exactly holds, the E-scheme reduces into the lattice scheme. 

Integrating Eq. (fT2l i. we obtain the well-known two-loop asymptotic scaling relation, 



R(gL/E) = «(Sl/e)A l/e = (b 0g t /E y bl/(2 ^exp 



-1 



2b g L/E 



(15) 
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Figure 6: Scaling at Nf = 6 from the N t dependence of the normalised critical temperature. Left: The bare lattice scheme results. The red symbol 
X shows T c /Al, and the blue □ symbols represent T c jh., obtained by using the parameters h = 0.03 and /3o = P£(Nf = 6, N, = 12) = 5.55 in 
Eq. Q6) , Right: The E-scheme results T c /Ae- 



with the lattice (E-scheme) Lambda-parameter, Al (Ae). In order to compensate the higher order loop effects, we 
consider the renormalisation group inspired improvement l6lll 



fl lmp (/?L/ E ) = A™ P E a(j3 L/E ) ee 



1 +h 



1 +h 



R 2 (J3l/e) 



(16) 



where /Jl/e = 10/(^l/e) 2 - The coupling /Jo will be specified to one of /? L C , and, the parameter h will be adjusted so 
as to minimise the scaling violation involved in /3 L C . Note that h — reproduces the standard asymptotic scaling law 
Eq. CTJ). 

The asymptotic scaling explained above is valid in the massless limit. In the following, we will use it to analyse 
results obtained at finite bare fermion mass ma = 0.02 by assuming that the shift of the (pseudo) critical coupling 
induced by a non-zero mass is smaller than other errors. This assumption should ultimately be tested in future studies 
by performing simulations with different masses and extrapolating to the chiral limit. 

In order to investigate the critical temperature T c in terms of our critical couplings /J L C or /J E C = 10/(gg) 2 , we 
substitute /J L C /E into the temperature definition Eq. ([TJi, and insert the scale Al/e: 



1 

N't 



A 



L/E 



: (a L /e a(0L/E)) • 



(17) 



The left-hand side is a given number, and Al/e a (P^/ E ) in the right-hand side is evaluated by using Eq. (TT3T > and 
our critical couplings fl£/E- Thus, Eq. (fTTT i allows us to convert the critical couplings into the (normalized) critical 
temperature r f /Al/e- When we adopt the improvement Eq. ( fl6l l. 7\ /Al/e is upgraded into r c /A™ E . 

4.2. Asymptotic scaling in terms of T c / Al/e and r c /A™ E 

The left panel of Fig.|6]shows T c /Kl as a function of the temporal extension N, at Nf = 6. As indicated by the red 
(x) symbols, T c /A^ at Nf = 6 varies within less than 10 percent in the whole range of N t , and stays almost constant 
between N, — 8 and 12. As shown in blue □ symbols in the left panel, N t dependences are suppressed by using the 
improved asymptotic scaling with h = 0.03 and /Jo = P^(N, = 12) = 5.55. We find that a small improvement h <k 1 
efficiently reduce the N, dependence of the normalised critical temperature. Furthermore in the E-scheme, the ratio 
T c /Ae without the improvement is almost N, independent as shown in the right panel of Fig. [6] These observations 
indicate that our /J L C at Nf = 6 has been determined near to the continuum limit, and the effects of scaling violations 
are in the controllable domain. 
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Figure 7: Scaling at Ne = 8 from the N, dependence of the normalised critical temperature. Left: The bare lattice scheme results. The red symbol 
X shows and the blue □ symbols represent r ( ./A™ p obtained by using the parameters h = 1.08 and /?o = P£(Nf = 8, N, = 12) = 4.34 

in Eq. 116) . Right: The E-scheme results. The red symbol X shows T c /Al, and the blue □ symbols represent r c /A™ p obtained by using the 
parameters h = 0.4 and j3 Q = P£(N f = 8, N t = 12) 5.94 in Eq. (TEi . 



Table 2: Summary of T c /Al and T c /A'™' 1 ' for various (Nf, N t ). The first (second) line at fixed (Nf, N t ) shows the value of T c /Al (r c /A™ p ), and 
the last two columns provide the parameter h and j8o appeared in the improved asymptotic scaling Eq. )16t . 



N/\N, 


4 


6 


8 


12 


h 


Po 





18.11+0.65 


18.21+0.91 


16.56 + 0.71 










16.29 + 0.75 


17.81 + 1.02 


16.56 + 0.78 




0.05 


8.26 


4 


21.99+1.04 


19.98 + 0.95 


17.12 ±2.43 










16.56+1.44 


18.67+1.38 


17.12 ±3.41 




0.30 


6.15 


6 


25.41 ± 1.43 


25.33 ±1.43 


22.94 ± 1.29 


22.30 ±2.52 








21.66+1.64 


23.87+ 1.58 


22.21 ± 1.40 


22.30 ±2.66 


0.03 


5.55 


8 




50.05 + 0.87 


47.06 ± 3.28 


34.34 ± 1.91 










34.32+ 1.40 


42.67 ± 6.33 


34.34 ±3.90 


1.08 


4.34 
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Table 3: Summary of T c /A^ and r c /A™g for Nf = 6 and Nf = 8. The first (second) line at fixed (Nf,N t ) shows the value of T c /A^ (7" C /A™g), 
and the last two columns give the parameter h and /3o appeared in the improved asymptotic scaling Eq. Q6) , For Nf = 6, the improvement was not 
necessary. 



N f \N, 


4 


6 


8 


12 


h 


A> 


6 


74.22 ± 5.86 


75.47 + 8.17 


74.56 ± 9.08 


75.13 + 10.76 




8 




422.54 + 23.06 


422.61 ± 38.59 


316.03 + 20.06 










312.16 + 33.13 


393.58 ± 60.01 


316.03 ±31.52 


0.40 


4.34 



Now we shall investigate the asymptotic scaling for the Nf = 8 system. In the left panel of Fig. [7] we show the 
N, dependence of the normalised critical temperature in the lattice scheme. The red symbol x shows the ratio T c /Al, 
which shows about 30 percent variation between N t - 6 and 12. The blue □ symbols represent 7V/A™ P obtained by 
using the parameters h = 1.08 and /3o = P£(Nf = %,N t = 12) = 4.34 in Eq. ( TTBT l. Although we have performed a 
large improvement h ~ 0(1), the variation still remains at 20 percent level. It turns out that the use of larger h cannot 
reduce the variation. 

The right panel of Fig. [7] shows our E-scheme results. The red symbol x shows the ratio T c /Al. Now the N, 
dependence almost disappears between N t - 6 and 8. The blue □ symbols represent 7V/A™ P obtained by using the 
parameters h = 0.4 and fio = /3£(Nf = 8,N t = 12) = 5.94 in Eq. (fT6l. The N, dependences are slightly suppressed 
by the improvement, while it cannot make all three data being consistent. Again, a larger h cannot reduce the N, 
dependence. 

Comparing to Nj = 6 (or smaller as shown in followings) cases, the thermal scaling of y6 L c at Nj = 8 shows 
much larger and less controllable deviations from the two-loop asymptotic scaling. This means that our /? L C involves 
a physics which cannot be captured by using the two-loop asymptotic scaling. This might be natural, in view of the 
largish values of the coupling involved. In particular, an IRFP occurring at finite coupling requires corrections to two 
loop scaling: in other words, an almost exact two loop scaling for Nf = 8 would rule out an IRFP at finite coupling 

At the same time, since the first two coefficients are universal, we expect to recover asymptotic scaling at very 
large N,, and very small couplings, in the deep UV region 

We have also investigated the asymptotic scaling for Nf = and 4 systems. Table [2] provides the summary of 
T c /Al and J c /A™ p for all Nf = 0, 4, 6, and 8. The first and second lines in the summary table provide quench (Nf = 
0) results for N, = 4, 6, and 8. Here T c represents the critical temperature associated with the deconfinement transition 
rather than the chiral transition. We find about 10 percent variation of for T c /Al in the whole range N t e [4, 8]. The N, 
dependence can be reduced into less then 10 percent by using the improved asymptotic scaling with a small h = 0.05. 

The third and forth lines in the table Table[2]show our Nf — 4 results for T c /Al and T c /A™ p for N, = 4, 6, and 8. 
We find about 20 percent variations of T c /Al between N, = 4 and N, = 8 results. The improved asymptotic scaling 
Eq. ( TTBT l works well, and the variation reduces into 10 percent level in whole range of N, =4-8. 

In summary, T c /A computed using different schemes (A = Al or Ae) consistently shows an increase with Nf, 
confirming the findings of our early work fl, and indicating that Al/e as defined in this Section vanishes faster than 
T c upon approaching the critical number of flavour. Within the various uncertainties discussed here, this can be taken 
as a qualitative indication of a scale separation. 

5. Onset of the conformal window 

In this section, we investigate the vanishing of the chiral dynamics with increasing Nf by utilising our /? L C and 
discuss the emergence of the conformal phase. 

In the first Subsection l5.ll we investigate so-called Miransky-Yamawaki (Nf - g) phase diagram, and discuss how 
this diagram can be used to estimate the onset of the conformal window Af* from the vanishing thermal scalings of 
jS L c . In the second Subsection 15.21 we introduce the thermal critical coupling as a typical interaction strength at 
critical temperature T c by utilising our /? L C , following a recent work by Shuryak Jioll . Then we discuss the onset of 
the conformal window from the approach of gj to the zero T critical couplings at large Nf. In the final Subsection 
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we study the Nf dependence of the critical temperature by introducing a UV Nf independent reference scale, and 
estimate Nj from the vanishing critical temperature. 
Here, we summarise our estimate Ny. 



N} 



11.1 + 1.6 (from the vanishing thermal scaling of yS L c ) , 

12.5 ± 0.7 (from the approach of g^ to g Ct so and £irfp,4i) , (18) 
10.4 + 1.2 (from the vanishing of T c /M with a scale M given by u = 0.79 - 0.81) . 

Details will be explained in the following. 

5.1. Miransky-Yamawaki phase diagram with thermal step-scaling 

We first recall well-known step-scaling analysis, which is based on the uniqueness of the physical correlation 
length 

f = a08 L )| = a08 L ')l' ■ (19) 

Here £ is the dimension-less lattice correlation length. Then, the vanishing (smaller) step-scaling — /?l - Pl is 
interpreted as a signal for a (pre-)conformal fixed point. 

Similarly, the uniqueness of the critical temperature at fixed Nf 

T; 1 = N t fl06 L c ) = N' t fl(/3 L c ') , (20) 

with N, + N' t gives a step-scaling A/? L C = jS L c -jS L c '. Again, the vanishing (attenuating) A/3 L C would be the signal of the 
(pre-)conformality. 

We here estimate the number of flavour satisfying A/3^(N*j) = by extrapolating our /3 L C collection into the larger 

flavour region. To this end, we plot the (pseudo) critical lattice coupling g L c = J10/j3£ as a function of Nf in Fig. [8] 
which gives an extension of so-called Miransky-Yamawaki (MY) diagram to finite temperature cases. 

In Fig. [8] g£ = -^/10//? L C at Nf > 6 is plotted in the MY diagram. Here we shall consider two lines characterized 
by N t = 6 and 12. We find that the lines cross at (g£,Nj ) = (1.79 + 0.12, 1 1.1 + 1.6). The vanishing of the step scaling 
at this point indicates the emergence of the conformal phase. 

We underscore that at this stage this is merely a qualitative discussion : the precise shape of the lines at fixed N, is 
dictated by the beta-function, which is unknown. Hence, a linear extrapolation which only uses two values of N, has, 
at this stage, only the meaning to illustrate a viable strategy which we plan to further pursue in the future. This said, 
the agreement of the results found here with the estimates presented below is rather gratifying. 

Before closing this subsection, we comment on another interesting aspect of obtained critical lattice couplings 
/? L C . As indicated by the constant N, line in Fig. [8] the critical coupling is an increasing function of Nf for a fixed 
lattice temporal extension. This behavior is the direct consequence of enhanced fermion screening effects due to the 
larger number of fermion species. We have observed the thermal scalings for our /? L C , thereby, the enhancement of the 
screening effects gets to a physical significance. The study of the thermal scaling at large Nf has been pioneered by 



J. Kogut and his collaborators 16211 . and the present results extend this early work, putting it in connection with the 



physics of the conformal window. 

5.2. Conformal phase in terms of thermal critical coupling 

In this subsection, we investigate the vanishing nature of the critical temperature T c at large Nf by introducing the 
thermal critical coupling. This approach to the conformal dynamics is motivated by the recent work Ref. Jul. Let us 
consider the renormalisation flow by integrating the two-loop beta-function 

' dg L 



M(?, ret ) r f«i 



(g L c ) 2 (gfVh+b o)-"^ 
{gO 2 bx+b (£- f ) 2 



x exp 



2^ol(^' rf ) 2 i8D 2 ' 



(21) 
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Figure 8: (Pseudo) critical values of the lattice coupling g£ = ^JlQ//3£ for theories with Nf = 0, 4, 6, 8 and for several values of N, in the 

Miransky-Yamawaki phase diagram. We have picked up g L c at Nf = 6 and 8, and considered a "constant N," line with N t = 6, 12. If the system 
is still described by one parameter beta-function in this range of coupling, the IRFP could be located at the intersection of the fixed N, lines - or 
equivalently, in the region where the step scaling function vanishes. To demonstrate the procedure -as a preliminary example - we have considered 
the intersection of the Nt = 12 and Nf = 6 lines. 



Here, our critical coupling g L c = -y/lO//?^ is used as the UV initial condition in Eq. ( |2TI ). and the renormalisation 
flow ends up at a certain reference scale M(g^ et ) specified by giving the reference coupling g L ref . The limit t ^^ ef — > oo 
in Eq. (l2TT i reproduces the standard asymptotic scaling formula Eq. (T5[ upto the overall constant (to account for 
conventional A L ) and the higher order of g L c , R(g£,g^ ef — » oo) — » C(l + 0((g£) 2 ))R(g£) . 
For the UV initial condition in Eq. (|2TT >. we particularly pick up the lattice critical coupling 



8l = ^m.- 



1.100 ±0.004 (Nf 

1.275 ±0.040 (Nf 

1.342 ±0.032 (N f 

1.518 ±0.021 (N f 



as a best data set from our critical coupling collection: Table Q] 
Now we re-write the critical temperature in relation to M(g L ref ), 

T r 1 



0, N, 

4, N t 

6, N t 

8, N t 



1 



8) 
8) 
12) 
12) 



(22) 



M(g L ref ) 



N a(g°) x M(g , ref ) 



(23) 



The reference coupling g^ ef which makes the reference scale M(g L J ef ) be equivalent to r e in Eq. ( 1231 is interpreted as 
a typical interaction strength at T c , We call such a coupling the thermal critical coupling g^: T c /M(g^) = 1. 

We shall now investigate the onset of the conformal phase by utilising g T c . The red (□) symbol in Fig.|9]shows gj 
for various Nf, and the dashed red line represents the fit for the data by adopting the ansatz 



N f (g°) = A ■ log [B ■ (g< - g T V =0 ) + 1] . 



(24) 



Here, A and B are the fit parameters. We compare g T c with the zero temperature critical couplings (g c ,sn blue dashed 
line) estimated by the two-loop Schwinger Dyson equation |[50ll as well as the IRFP position (giRFp,4i, magenta dashed- 
dotted line) of the four-loop beta-function. For the latter, the central value is determined by using the MS scheme, 
and the error-bars represent a possible scheme dependence [49]. At the crossing points, the thermal and zero T chiral 
transition coincide, which indicate the vanishing critical temperature with the emergence of the conformal phase. 
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g 

Figure 9: The thermal critical coupling (red □) and the fit for them (dashed red line, with the ansatz Eq. (24)) and the values of the zero temperature 
couplings in the conformal phase from different estimates, see text for details. At the critical number of flavour the thermal critical coupling should 
equal the critical coupling associated with the IRFP. The procedure is motivated by a recent study by Shuryak in Ref. IToh . 



Thus the candidate for the IR fixed point is estimated to be (g*,Nj) = (2.79, 13.2) or (2.51, 1 1.8) from the joint of the 
gj -gc.sD or the g T c -£irfp,4i, respectively. Combining both results, we get Nj ~ 12.5+0.7. Besides the conformality, 
another interesting property of g T c is its increasing trend at larger Nf, which indicates a realisation of more strongly 
interacting non-Abelian plasma at larger Nf as discussed in Ref. ifTofl . 

5.3. Conformal phase in terms of vanishing critical temperature 

We investigate the onset of the conformal window in terms of the vanishing critical temperature T c with increasing 
Nf. The pioneering study on the Nf dependence of T c is Ref. Jit], where the functional renormalisation group method 
was adopted. Also, T c {Nf) has been studied in the bottom-up approach based on the gauge-gravity duality Q]. The 
first attempt along the course of this direction in the lattice framework is found in our previous work Jgt], where we set 
a reference scale M by specifying a certain value for g^ e{ in Nf independent way in Eq. d2T1 i. and investigated the Nf 
dependence of the ratio T c /M. In this subsection, we elaborate on the way to introduce a relevant reference scale M by 
utilising zero T plaquette data as well as the thermal critical coupling g^. Before going to details, we first explain the 
basic idea by quoting the the FRG study by Braun and Gies. They used the r lepton mass m T = 1 .777 (GeV) as a Nf 
independent UV reference scale for theories with any number of flavours ifj. More specifically, the initial condition 
of the renormalization flow has been specified via the strong coupling constant in Nf independent ways: 

a s (p = m T ) = 0.322 for *N f . (25) 

Starting from the common initial condition Eq. d25l l, a Nf dependent screening strength associated with a Nf fermion 
multiplicity gives rise to a different pattern of the FRG renormalization flow at a chiral phase transition scale p ~ 
Aqcd ^ m T , giving a Nf dependent nature for a critical temperature T c (Nf). Then, the vanishing T c (Nf) with 
increasing Nf indicates a loss of the typical scale Aqcd, *•<?• the emergence of the conformality. Most importantly, 
the Nf dependence of T c as well as its novel non-analytic behavior in the pre-conformal region becomes free from the 
choice of the reference scale |Ql by taking an Nf independent UV reference scale much larger than T c . 

We here try to perform an analogous study to the FRG in the lattice framework within the available data. Ideally 
speaking, we would like to set the scale by measuring on the lattice some physical quantity insensitive to IR dynamics 
- for instance fixing the value of a s in the V-scheme to some appropriate value, as done in the computation of ro or 
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Figure 10: The /?l dependences of the tadpole factor uq at zero temperature (12 lattice volume). At each Nf, the dashed line represents the fit 
for data with the ansatz «o = 1 — A/(l + B ■ Bj). We consider a constant uq (e.g. uq = 0.8 in figure), and read off the corresponding lattice bare 
couplings which are used to define the scale M at each theory with Nf flavours. 



variations thereof, following Ref. 1 59] and related applications. These large scale simulations are now starting 165H . 
Here, we approximate the renormalisation flow by Eq. d2"TT i and d23l ). 

In order to determine the reference coupling g/ e which appears as in Eq. (I2H . we utilise our plaquette results 
(tad-pole factor uq) data shown in the left panel of Fig. Q~T] Let us consider a constant uq without Nf dependences, 
for instance uq = 0.8 in figure, and read-off the corresponding bare lattice couplings at each Nf. The obtained 
gh(Nf) = -y/10 1 'PiXN '/) is used as a reference coupling g^f in Eq. (1231) . 

Some remarks on the aforementioned scale setting are in order: First, we recall the scale setting procedure in the 
potential scheme, where the measured normalised force r 2 F(r) is proportional to the renormalised coupling g, and 
the specification g 2 oc r^F(rx) = 3 X sets a scale r x l . In short, we use our uq (or equivalently plaquettes) as g, and 
mo = X is regarded as the analog of the potential scheme scale setting. Second, in the leading order of the perturbative 
expansion, the renormalised coupling is Nf independent, and proportional to the Wilson loop lr63ll a property that we 
have already exploited in the E-scheme calculation. Hence the use of Nf independent «o approximately gives a Nf 
independent scale setting, similarly to the FRG scale setting method Eq. (|25T >. And third, such an Nf independent scale 
setting can be performed at a sufficiently UV regime T c (Nf) <k M(g^ e{ ) by adjusting the value of uq = X to satisfy 
the condition g/ 6 <k gjCNf). Here, the thermal critical coupling g T c has already been evaluated in the previous 
subsection. 

We note that the critical lattice coupling /? L C itself is NOT the coupling at the scale of T c = [N t ■ a(gf)]~ l , but the 
coupling at the lattice cut-off a~ l (g£), which is N, » 1 times larger than T c . Then, the scale hierarchy T c (Nf) <K 
a_1 (§L (Nf)) allows us to consider the reference scale much larger than critical temperature but smaller than the lattice 
cutoff T c (Nf) «: M(gl ef ) < a~\g£(Nf)) by imposing the condition g£(Nf) < g L ref «: gj(N f ). We find that u ~ 0.8 
meets this requirement. In summary, the use of g L given by uq ~ 0.8 would be analogous to the FRG scale setting 
method Eq. d25l l, and suitable for studying the vanishing chiral dynamics by utilising T c /M(g£ et ). For the UV initial 
condition in Eq. (fST), we use the values listed in Eq. d22l >. 

The middle panel of Fig. QT] displays the N f dependence of r c /M(g L ref ) defined by Eq. d23j for u = 0.79, 0.80, 
and 0.81. All T c /M(g^ ef ) data for Nf > 4 turn out to be smaller than 0.4, hence the scale M(g^" ef ) is at least two times 
larger than the chiral transition scale. Fitting the data points for T c /M(g£ et ) at Nf > 4 by using the FRG motivated 
ansatz, 

T c = K\N} - N f \ { - 2b ° /bi)(N r } , (26) 
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Figure 11: Left: The Nf dependence of T c /M where M is determined to be a UV scale corresponding to itn = 0.79 (red box), 0.80 (blue O). 
and 0.81 (magenta triangle) at each theory with Nf. The dashed lines represent fits for data by assuming the expression Eq. i26\ . Right: The uq 
dependence of N^. The three data in the left side are determined within the condition M(g L ref ) < a~'(g L c ), while for the others M(g^ f ) exceeds 
the lattice cutoff. This more robust procedure confirms our early results, and should be ultimately confirmed by use a rigorous lattice scale setting 
which is in progress loHl . 



where &o,i has been defined in Eq. ([TBI and (1141) . the lower edge of the conformal window is estimated as: Nj = 
10.48 + 1.01 (m = 0.79), N* = 10.34 + 0.88 (u = 0.80), N* = 10.23 + 0.80 (u = 0.81). The error-bars involve both 
fit errors and statistical errors of data. Combining them, we find N f ~ 10.4 + 1.2. 

As shown in the right panel of Fig. QT| Nj is relatively stable within the range 0.79 < uq < 0.94 though uq > 0.81 
gives a too large reference scale M(gi) > a^ 1 (g^) for some Nf, and then we should expect considerable discritization 
errors in addition to those we have presented in the figure. A small uq < 0.7 leads to M(g^ e{ ) ~ T c or smaller, then the 
decreasing nature of T c /M(g™ f ) is lost. In such a case, the reference scale M(g^ ef ) as well as T c is expected to include 
the chiral dynamics, and the vanishing of the chiral dynamics becomes no longer visible in the ratio T c /M(g^ et ). The 
decreasing T c /M(g™ f ) with UV M(g^ et ) and the loss of the decrease with non-UV M(g L ref ) are qualitatively consistent 
to the FRG observations |@]. 



6. Summary 

We have investigated the chiral phase transition (crossover) at finite temperature T in colour SU(3) QCD-like 
theories with various number of flavours Nf — (quench), 4, 6, and 8 by using lattice Monte-Carlo simulations. We 
have used a single bare fermion mass ma = 0.02. For all number of flavours, we have used the Asqtad action with 
a one-loop Symanzik and tadpole improved gauge action. The main focus in this paper is to investigate the chiral 
crossover at finite T as a function of Nf, and discuss the possible implication for the (pre-)conformal dynamics at 
large Nf. In Eq. dl8t . we provide the summary for the number of flavour Nl where the conformal window would 
emerge. 

The observables in our simulations were the chiral condensate, the Polyakov loop, and their susceptibilities for 
various lattice couplings fl^, lattice sizes, and the number of flavours Nf. And then, we have collected the (pseudo) 
critical lattice coupling as a function of (Nf, N t ). Table[T]provides the summary for /? L C . Our /? L C are consistent with 
enhanced fermionic screening at larger Nf. The use of several Af, allows us to perform the asymptotic scaling analyses 
for the thermal scaling of y6 L c . Further utilising y6 L c , we have discussed a possible implication for the (pre-)conformal 
dynamics at large Nf. 

We have estimated the Nl from the vanishing thermal scaling by extrapolating our critical couplings g L ° to larger 
Nf direction. This gives Nj ■ ~ 11.1 + 1.6. We have extracted a typical interaction strength gj at the scale of critical 
temperature T c by utilising our g L c and the two-loop beta function, and compared g T c to the zero temperature critical 
couplings (g c ,SD) estimated by the two-loop Schwinger Dyson equation 1 50] as well as the IRFP position (§irfp,4i) 



of the four-loop beta-function [49]. The coincidence between g T c and g e ,SD or £irfp,4i indicates the vanishing critical 
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temperature with the emergence of the conformal phase. Based on this reasoning, we have estimated the onset of the 
conformal window as ~ 12.5 ± 0.7. We have also confirmed the increasing of gj at larger Nf which has been 
discussed in Ref. fioll and indicates more strongly interacting non-Abelian plasma at larger Nf. 

Further, we have examined the Nf dependence of T c /M by introducing a UV Nf independent reference scale M 
which is determined by utilising the tadpole factor uq in analogous ways to the potential scheme scale setting. Then, 
T c /M turns out to be a decreasing function of Nf consistently to the FRG observations Jg], and the vanishing T c /M 
indicates the emergence of the conformal window around N^ ~ 10.4 ± 1.2. 

As a future perspective, we plan perform more rigorous scale settings, by exploiting state-of-art measurements of 
lattice potential. It is also mandatory to investigate the chiral limit and the thermodynamic limit at large Nf. This, 
together with a more extended set of flavour numbers, will allow a quantitative analysis of the critical behaviour in the 
vicinity of the conformal IR fixed point. We expect that our thermodynamic lattice study for the large Nf non-Abelian 
gauge theory plays an important role as a new connection between the lattice and the Gauge/Gravity duality l64ll7ll. 
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Appendix A. Summary table for simulation results 

In this appendix, we summarise our results and the parameters used in the simulations and analyses. 
We have used an improved version of the staggered action, the Asqtad action (52], with a one-loop Symanzik 
J53, 54] and tadpole |55j] improved gauge action. All our simulations used the same bare fermion mass ma — 0.02. 



In the following tables we will quote the lattice bare couplings /3l = lO/g^, tadpole factors uq, step lengths for a 
single trajectory At = 5t x 20, the number of total trajectories n tra j, the number of thermalised trajectories n ave which 
have been used to evaluate ensemble averages, bin-sizes for jackknife error estimates st,in> ensemble averages for the 
chiral condensate (2a 3 (ifnfr) with the definition Eq. (0)) the Polyakov loop (Re L or \L\ with the definition Eq. ([Toil), 
and the chiral susceptibility (Xa-, defined in Eq. (|7]i) and/or the susceptibility (x\l\) associated with the absolute value 
of Polyakov loop. 

Here are several technical comments for the symbols: In the second column in each table, some values for the 
tadpole factors uo have been estimated by using the fit with the ansatz uq = 1 .0 - A/(l .0 + B j3^), and are shown inside 
of the parentheses. Some of the total trajectories in the third column have a symbol +, for which the configurations 
obtained in the other simulations has been utilised as the inputs. The Monte-Carlo step At in the fourth column 
has been adjusted to give about 75 - 80 percent Metropolis acceptances in the pre-thermals domain. In the fifth 
column, some entries have a symbol *, which indicates that additional simulations would be preferable to confirm the 
thermalisation though the trajectories seem to have reached a stable domain. In the seventh column, a pre-factor 2 
appears in front of a 3 {ifitfr) following the difference between the the MILC code convention and the standard definition 
Eq. ©. 
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Table A. 4: Summary table for Nf = with the use of the lattice volume 16 3 X 4. 





«0 


Itraj 


At 


**ave 


■^bin 




\L\ 


X\L\ 


6.80 


(0.864) 


4000 


0.60 


2000 


200 


0.575251(110) 


0.024(1) 


0.6 + 0.04 


7.00 


(0.868) 


4000 


0.60 


2000 


200 


0.559423(113) 


0.032(2) 


1.10 + 0.11 


7.10 


(0.871) 


4000 


0.60 


2000 


400 


0.552847(121) 


0.038(7) 


2.29+1.10 


7.15 


(0.8715) 


3000 


0.60 


2000 


250 


0.549450(102) 


0.050(3) 


2.27 + 0.28 


7.20 


(0.872) 


4000 


0.60 


2000 


400 


0.546144(63) 


0.052(7) 


2.99 + 0.41 


7.25 


(0.8735) 


3000 


0.60 


1500 


250 


0.543789(46) 


0.051(4) 


2.06 + 0.32 


7.30 


(0.875) 


4000 


0.60 


1500 


250 


0.541361(119) 


0.057(10) 


4.75 + 0.65 


7.35 


(0.876) 


12000 


0.70 


6000 


1200 


0.538570(127) 


0.273(29) 


29.34+ 10.51 


7.40 


(0.877) 


12000 


0.68 


6000 


1200 


0.535896(83) 


0.423(7) 


5.43 + 0.90 


7.45 


0.878 


8000 


0.68 


4000 


500 


0.531813(120) 


0.395(10) 


9.59+1.83 


7.50 


(0.879) 


4000 


0.60 


2000 


400 


0.529767(146) 


0.440(14) 


7.30 + 0.66 


8.00 


(0.889) 


1000 


0.60 


500 


100 


0.516562(107) 


0.702(7) 


2.51+0.69 



Table A. 5: Summary table for Nf = with the use of the lattice volume 16 3 X 6. In order to find the suitable At, we performed some preceded 
simulations with the order of 5000 trajectories, and their output configurations have been used as inputs to get the results in this table. 





M 


«traj 


At 


^ave 


■^bin 


2a 3 (ifnf/) 


\L\ 


X\L\ 


7.50 


(0.879) 


4000+ 


0.68 


600 


100 


0.532814(60) 


0.0278(018) 


0.77 + 0.08 


7.60 


(0.881) 


8000+ 


0.68 


1000 


200 


0.529409(42) 


0.0372(018) 


1.72 + 0.32 


7.70 


(0.883) 


12000+ 


0.68 


2000 


400 


0.526335(44) 


0.0652(086) 


5.34 + 0.75 


7.80 


(0.885) 


18000+ 


0.68 


2000 


500 


0.523805(11) 


0.0656(107) 


4.44+1.41 


7.90 


0.886994(54) 


10000+ 


0.68 


2000 


400 


0.521344(35) 


0.1950(157) 


8.20 + 3.19 


8.00 


0.889 


6000+ 


0.68 


900 


100 


0.519265(44) 


0.2538(079) 


5.08+0.51 


8.10 


0.890728(45) 


6000+ 


0.68 


1000 


250 


0.517268(47) 


0.2838(085) 


4.48+0.51 


8.20 


0.892579(35) 


6000+ 


0.68 


1000 


125 


0.515598(21) 


0.3255(042) 


3.43 + 0.28 


8.30 


0.894318(46) 


4000+ 


0.64 


600 


100 


0.514019(24) 


0.3514(061) 


3.93 + 0.42 


8.40 


(0.896) 


4000+ 


0.64 


600 


120 


0.512534(31) 


0.3806(049) 


3.43 + 0.55 



Table A. 6: Summary table for Nf = with the use of the lattice volume 24 3 X 8. 





u 


«traj 


At 


'^ave 


■^bin 


2a 3 (iAiA) 


\U 


X\L\ 


7.80 


(0.885) 


4800 


0.52 


1000 


200 


0.523836(26) 


0.0107(5) 


0.37 + 0.02 


8.00 


(0.889) 


3100 


0.56 


2000 


400 


0.519532(22) 


0.0139(10) 


0.68 + 0.03 


8.10 


0.890728(45) 


12000 


0.54 


1000 


250 


0.517575(22) 


0.0175(21) 


1.00 + 0.09 


8.15 


(0.892) 


9600 


0.56 


4000 


800 


0.516582(15) 


0.0566(64) 


4.31 + 1.06 


8.20 


0.892579(35) 


8800 


0.52 


2000 


400 


0.515820(23) 


0.0317(32) 


2.19 + 0.45 


8.25 


0.893 


7200 


0.54 


3000 


600 


0.514974(12) 


0.0770(71) 


4.89+1.24 


8.30 


0.894318(46) 


8700 


0.56 


3000 


600 


0.514153(9) 


0.1158(43) 


3.32 + 0.94 


8.40 


(0.896) 


4800 


0.56 


1500 


300 


0.512647(12) 


0.1623(36) 


2.23 + 0.47 


8.60 


(0.899) 


4000 


0.56 


1500 


300 


0.509860(4) 


0.1807(45) 


2.63 + 0.79 


8.80 


(0.900) 


4000 


0.56 


1500 


300 


0.505961(5) 


0.2160(39) 


3.53 + 0.42 


9.00 


(0.900) 


4000 


0.56 


1500 


300 


0.501632(5) 


0.2587(33) 


2.41+0.56 
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Table A. 7: Summary table for Nf = 4 with the use of the lattice volume 16 3 X 4. 



A. 


u 


"traj 


At 


^ave 




2<rW> 


X<T 


Re L 


4.80 


(0.774) 


3000 


0.12 


1000 


250 


0.5614(15) 


0.356(17) 


0.0053(21) 


5.10 


(0.796) 


4000 


0.12 


1500 


250 


0.4733(13) 


0.552(17) 


0.0141(16) 


5.30 


(0.809) 


4000+ 


0.14 


1500 


250 


0.4001(13) 


0.688(22) 


0.0256(30) 


5.40 


(0.815) 


6000 


0.14 


1500 


300 


0.3485(26) 


0.934(49) 


0.0406(19) 


5.45 


(0.818) 


5000+ 


0.14 


3000 


500 


0.3132(26) 


1.090(77) 


0.0579(14) 


5.50 


(0.821) 


4000+ 


0.14 


1500 


250 


0.2893(11) 


0.993(24) 


0.0657(15) 


5.55 


(0.824) 


5000+ 


0.16 


2000 


250 


0.2443(16) 


1.281(48) 


0.0917(23) 


5.60 


(0.826) 


7000 


0.16 


2000 


500 


0.1814(25) 


1.746(103) 


0.1356(35) 


5.65 


(0.829) 


5000+ 


0.18 


2000 


400 


0.1117(26) 


1.934(97) 


0.1953(41) 


5.70 


(0.831) 


4000 


0.20 


1500 


250 


0.0859(14) 


1.747(15) 


0.2231(27) 


5.80 


0.836159(98) 


4000+ 


0.24 


1500 


250 


0.0610(7) 


1.425(8) 


0.2590(23) 


6.00 


0.846224(81) 


3000 


0.28 


1500 


150 


0.0391(1) 


0.964(2) 


0.3519(28) 



Table A. 8: Summary table for Nf = 4 with the use of the lattice volume 16 3 X 6. 



A- 


M 


w traj 


At 


^ave 




2a\W) 


Xa 


Re L 


5.70 


(0.831) 


3500 


0.12 


1000 


200 


0.2311(8) 


0.948(8) 


0.0025(15) 


5.80 


0.836159(98) 


4500 


0.16 


1000 


200 


0.1772(13) 


1.102(29) 


0.0125(10) 


5.85 


(0.8385) 


6000 


0.16 


1000 


200 


0.1617(9) 


1.057(10) 


0.0094(19) 


5.90 


(0.841) 


5500 


0.16 


3000 


750 


0.1338(28) 


1.270(94) 


0.0180(32) 


5.95 


(0.8436) 


6000 


0.16 


2000 


500 


0.1000(14) 


1.311(30) 


0.0394(18) 


6.00 


0.846224(81) 


8000 


0.18 


4000 


800 


0.0831(25) 


1.499(98) 


0.0463(33) 


6.10 


0.850711(61) 


2500 


0.20 


1000 


200 


0.0503(4) 


1.152(12) 


0.0854(26) 


6.20 


0.854878(58) 


1600 


0.20 


600 


100 


0.0413(3) 


1.003(7) 


0.1081(20) 


6.30 


(0.859) 


1500 


0.24 


500 


85 


0.0366(2) 


0.902(5) 


0.1198(21) 


6.40 


0.862358(65) 


1200 


0.24 


400 


80 


0.0330(2) 


0.816(33) 


0.1464(28) 
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Table A. 9: Summary table for Nt =4 with the use of lattice volume 24 3 X 8. For the simulations with /?l = 6.5 and 6.8, the Monte-Carlo step Ar 
must be set to be smaller than the presented values in the early stage of the molecular dynamics evolutions to avoid low Metropolis acceptances. 





«0 


"traj 


Ar 


^ave 


*bin 


2fl J (i/r(A) 


Xo- 


Re L 


5.60 


(0.826) 


1500+ 


0.12 


400 


100 


0.2829(8) 


0.861(23) 


-0.0004(23) 


5.80 


0.836159(98) 


2100 


0.14 


600 


75 


0.1823(3) 


0.965(19) 


0.0009(8) 


5.90 


0.841 


1500+ 


0.14 


1000 


200 


0.1453(5) 


1.017(12) 


0.0014(5) 


5.95 


(0.8436) 


4000+ 


0.14 


1600 


400 


0.1268(10) 


1.176(26) 


0.0015(6) 


6.00 


0.846224(81) 


7000 


0.16 


2500 


625 


0.1098(10) 


1.113(24) 


0.0031(3) 


6.05 


(0.848) 


4800+ 


0.16 


3500* 


875 


0.0973(11) 


1.107(9) 


0.0046(5) 


6.10 


0.850711(61) 


5400 


0.16 


2000 


500 


0.0828(7) 


1.066(23) 


0.0068(4) 


6.15 


(0.853) 


2000+ 


0.16 


800 


100 


0.0690(3) 


1.046(10) 


0.0122(9) 


6.20 


0.854878(58) 


2500 


0.18 


800 


114 


0.0584(3) 


1.046(9) 


0.0165(4) 


6.25 


(0.8567) 


4000+ 


0.18 


2000 


400 


0.0518(3) 


1.035(13) 


0.0209(13) 


6.30 


(0.859) 


2540 


0.22 


1000 


200 


0.0435(3) 


1.006(7) 


0.0280(13) 


6.35 


(0.8608) 


1590+ 


0.22 


1200 


240 


0.0401(4) 


0.962(5) 


0.0346(10) 


6.40 


0.862358(65) 


3500 


0.24 


1500 


375 


0.0384(3) 


0.917(6) 


0.0378(8) 


6.45 


(0.8647) 


1590+ 


0.26 


800 


160 


0.0358(1) 


0.871(2) 


0.0425(13) 


6.50 


(0.867) 


2350 


0.24 


700 


100 


0.0340(1) 


0.830(2) 


0.0501(12) 


6.80 


(0.880) 


3300 


0.40 


1000 


200 


0.0284(0) 


0.707(1) 


0.0817(14) 



Table A. 10: Summary table for Nf = 6 with the use of the lattice volume 16 3 X 4. 





MO 


n traj 


At 


'^ave 




2a j <</^> 


Xo- 


Re L 


4.00 


(0.759875) 


2000 


0.10 


1000 


300 


0.5829(1) 


0.304(5) 


0.0011(18) 


4.30 


0.776360(228) 


3000 


0.10 


1500 


300 


0.4957(21) 


0.585(48) 


0.0151(36) 


4.40 


0.783740(347) 


3000 


0.10 


1500 


300 


0.4565(11) 


0.615(7) 


0.0204(7) 


4.50 


0.788558(123) 


4000 


0.10 


1000 


300 


0.4077(24) 


0.778(45) 


0.0349(48) 


4.60 


0.795206(97) 


5000 


0.10 


2000 


300 


0.3476(15) 


0.982(27) 


0.0557(14) 


4.65 


(0.798) 


4000+ 


0.16 


2000 


1000 


0.2667(52) 


2.609(451) 


0.0940(37) 


4.70 


0.800839(277) 


4000 


0.24 


1000 


300 


0.1244(22) 


2.206(38) 


0.1797(23) 


4.80 


0.805839(364) 


4000 


0.10 


1000 


300 


0.0822(12) 


1.831(15) 


0.2195(14) 


5.00 


0.817551(86) 


2000 


0.10 


1000 


300 


0.0521(6) 


1.271(8) 


0.2850(26) 


5.20 


0.828421(354) 


2000 


0.10 


1000 


300 


0.0418(3) 


1.035(8) 


0.3300(57) 


5.50 


0.841873(95) 


2000 


0.10 


1000 


300 


0.0325(0) 


0.810(1) 


0.4070(41) 
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Table A. 1 1 : Summary table for Nf = 6 with the use of the lattice volume 1 6 3 X 6. 





MO 


n traj 


At 


'^ave 


■^bin 




Xo- 


Re L 


3.60 


0. /3403 1(104) 


1 300 


0.16 


1 AAA 

1000 


OCA 

250 


0.6622(1) 


A 1 A / 1 \ 
0.14(1) 


A AAQ A t *Z\ 

—0.0034(6) 


4.00 


(0.759875) 


1500 


0.16 


1000 


250 


0.5923(5) 


0.30(1) 


-0.0041(4) 


4.20 


(0.773118) 


3000 


0.16 


1000 


250 


0.5422(7) 


0.55(14) 


-0.0021(23) 


4.40 


0.783740(347) 


3000 


0.16 


1000 


250 


0.4808(0) 


0.72(16) 


-0.0050(4) 


4.50 


0.788558(123) 


3000 


0.16 


1000 


250 


0.4363(2) 


1.12(26) 


-0.0032(11) 


4.60 


0.795206(97) 


3000 


0.18 


1000 


250 


0.3947(20) 


0.96(17) 


-0.0008(4) 


4.70 


0.800839(277) 


3000 


0.20 


1000 


250 


0.3418(10) 


1.12(17) 


-0.0007(7) 


4.80 


0.805839(364) 


3000 


0.24 


1000 


250 


0.2864(42) 


1.43(18) 


0.0028(3) 


4.90 


0.811809(354) 


3000 


0.20 


1000 


250 


0.2231(34) 


1.64(15) 


0.0058(14) 


4.95 


(0.81466) 


4000 


0.24 


1000 


250 


0.1935(17) 


1.37(5) 


0.0091(1) 


5.00 


0.817551(86) 


7000 


0.20 


2000 


500 


0.1644(28) 


1.66(8) 


0.0137(10) 


5.05 


(0.8196) 


6000+ 


0.24 


1000 


250 


0.1072(63) 


1.98(19) 


0.0388(32) 


5.10 


(0.821629) 


3000 


0.20 


1000 


250 


0.0769(21) 


1.51(10) 


0.0566(17) 


5.20 


0.828421(354) 


3000 


0.24 


1000 


250 


0.0581(2) 


1.38(1) 


0.0710(2) 


5.30 


0.832865(89) 


3000 


0.20 


1000 


250 


0.0489(0) 


1.13(6) 


0.0850(18) 


5.50 


0.841873(95) 


3000 


0.28 


1000 


250 


0.0393(0) 


0.97(0) 


0.1138(3) 



Table A. 12: Summary table for Nf = 6 with the use of the lattice volume 24 3 X 8. 





MO 


«traj 


At 


^ave 


^bin 


2a\W) 


Xa 


Re L 


4.60 


0.795206(97) 


1000 


0.20 


500 


100 


0.3985(4) 


0.68(4) 


-0.0014(11) 


4.80 


0.805839(364) 


1000 


0.20 


500 


100 


0.2929(15) 


1.08(5) 


-0.0008(0) 


4.90 


0.811809(354) 


1210 


0.20 


500 


50 


0.2335(16) 


1.22(12) 


0.0001(3) 


5.00 


0.817551(86) 


3500+ 


0.24 


1000 


200 


0.1840(7) 


1.31(2) 


-0.0001(4) 


5.10 


(0.821629) 


3000 


0.24 


1500 


300 


0.1357(8) 


1.30(3) 


0.0022(4) 


5.15 


(0.825) 


5000+ 


0.28 


2000 


400 


0.1164(10) 


1.37(6) 


0.0037(3) 


5.20 


0.828421(354) 


5000 


0.28 


2500* 


500 


0.0959(18) 


1.66(14) 


0.0067(5) 


5.25 


(0.8306) 


5000+ 


0.32 


3000 


500 


0.0732(9) 


1.50(7) 


0.0126(7) 


5.30 


0.832865(89) 


4000 


0.32 


1500 


300 


0.0606(5) 


1.36(1) 


0.0177(4) 


5.50 


0.841873(95) 


1500 


0.32 


500 


125 


0.0441(1) 


1.07(0) 


0.0301(7) 
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Table A. 1 3 : Summary table for Nf = 6 with the use of the lattice volume 24 3 X 1 2. 



A. 


«0 


Itraj 


At 


^ave 


■^bin 


2a J <i/^> 


X<T 


Re L 


4.70 


0.800839(277) 


1000 


0.20 


400 


100 


0.3493(6) 


0.753(20) 


-0.000188(377) 


4.80 


0.805839(364) 


1000 


0.20 


400 


100 


0.2961(4) 


0.927(35) 


-0.000096(510) 


4.90 


0.811809(354) 


1000 


0.20 


400 


100 


0.2330(7) 


1.200(70) 


-0.000344(692) 


5.00 


0.817551(86) 


1500 


0.20 


1500 


500 


0.1838(9) 


1.356(50) 


-0.000616(710) 


5.10 


(0.821629) 


1500 


0.20 


900 


300 


0.1420(12) 


1.269(63) 


-0.000297(360) 


5.20 


0.828421(354) 


1500 


0.20 


600 


200 


0.1065(2) 


1.172(23) 


0.000010(144) 


5.30 


0.832865(89) 


1500 


0.24 


900 


225 


0.0838(4) 


1.140(14) 


-0.000158(387) 


5.35 


(0.835280) 


1500+ 


0.28 


800 


200 


0.0748(4) 


1.138(16) 


0.000147(336) 


5.40 


(0.837659) 


2600+ 


0.30 


1600 


400 


0.0677(5) 


1.120(14) 


0.000185(79) 


5.45 


(0.839960) 


3100+ 


0.30 


2000 


400 


0.0622(2) 


1.060(7) 


0.000071(234) 


5.475 


(0.841081) 


2450+ 


0.32 


1500 


300 


0.0580(3) 


1.063(7) 


0.000202(188) 


5.50 


0.841873(95) 


4100+ 


0.32 


1500 


375 


0.0570(5) 


1.064(10) 


0.001716(248) 


5.525 


(0.843350) 


2600+ 


0.34 


1500 


300 


0.0536(2) 


1.041(7) 


0.000205(467) 


5.55 


(0.844445) 


2600+ 


0.36 


1500 


300 


0.0508(2) 


1.031(3) 


0.000545(238) 


5.575 


(0.845519) 


2450+ 


0.38 


1500 


300 


0.0483(4) 


1.073(17) 


0.000956(218) 


5.60 


(0.846582) 


4100+ 


0.34 


1000* 


250 


0.0456(4) 


1.040(12) 


0.001615(151) 


5.70 


(0.850883) 


900+ 


0.36 


400 


100 


0.0409(1) 


0.960(2) 


0.002867(470) 



Table A. 14: Summary table for Nf = 8 with the use of the lattice volume 24 3 X 8. The tadpole factors ug have been computed in I2H1 . 



A- 


M() 


Itraj 


At 


^ave 


*bin 


2a 3 <i/^> 




Re L 


4.00 


(0.793260) 


4500 


0.12 


1000 


200 


0.1584(8) 


2.72(3) 


0.00006(113) 


4.10 


0.79825 


3800 


0.12 


1200 


300 


0.1387(9) 


2.52(4) 


-0.00074(82) 


4.15 


0.80175 


2400+ 


0.14 


1000 


200 


0.1319(3) 


2.39(1) 


0.00005(92) 


4.20 


0.8053 


6500 


0.14 


1500 


375 


0.1232(7) 


2.30(5) 


-0.00094(85) 


4.25 


0.8086 


8700+ 


0.14 


8000* 


2000 


0.1069(18) 


2.06(15) 


0.00584(42) 


4.30 


(0.8116) 


2570 


0.16 


1000 


200 


0.0767(4) 


1.63(1) 


0.01282(50) 


4.40 


0.8192 


2570 


0.18 


1000 


200 


0.0639(2) 


1.48(0) 


0.01801(52) 
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